in Proceedings of the ECAI Workshop on Spatial and Temporal Reasoning, pp. 135-139, Valencia, Spain, 2004: 



A TCSP-like decidable constraint language generalising 
existing cardinal direction relations 

Amar Isli 
FB Informatik, Universitat Hamburg 
am99i@ yahoo.com 



WORK EXACTLY AS REJECTED AT THE MAIN ECAI 2 2004 3 
Abstract. We define a quantitative constraint language subsum- 
ing two calculi well-known in QSR 4 : Frank's cone-shaped and 
projection-based calculi of cardinal direction relations. The lan- 
guage is based on convex constraints of the form (a, j3)(x, y), with 
a, (3 £ [0,27r) and (fi— a) G [0, 7r): the meaning of such a constraint 
is that point x belongs to the (convex) cone-shaped area rooted at y, 
and bounded by angles a and (5. The general form of a constraint 
is a disjunction of the form [(ax, pi) V • • • V (a„, (3 n )](x,y), with 
(ai,f3i)(x, y), i = 1 . . .n, being a convex constraint as described 
above: the meaning of such a general constraint is that, for some 
i = l...n, (oi, /3i)(x,y) holds. A conjunction of such general 
constraints is a rCSPTike CSP, which we will refer to as an SCSP 
(Spatial Constraint Satisfaction Problem). We describe how to com- 
pute converse, intersection and composition of SCSP constraints, 
allowing thus to achieve path consistency for an SCSP. We show 
how to translate a convex constraint into a conjunction of linear 
inequalities on variables consisting of the arguments' coordinates. 
Our approach to effectively solving a general SCSP is then to 
adopt a solution search algorithm using (1) path consistency as the 
filtering method during the search, and (2) the Simplex algorithm, 
guaranteeing completeness, at the leaves of the search tree. 

Keywords: Constraint Satisfaction, Spatial reasoning, Geomet- 
ric Reasoning, Knowledge Representation, Qualitative Reasoning, 
Quantitative Reasoning 

1 Introduction 

Concilialing qualitative reasoning and quantitative reasoning in KR&R systems: 
a way lo systems representalionally more flexible, cognilively more plausible, and. eompulalionally. wilh the advantage of 
having the choice between a purely-quantitative and a qualilulive-compulalions-lirsl behaviours 

Knowledge representation (KR) systems allowing for the represen- 
tation of both qualitative knowledge and quantitative knowledge are 
more than needed by modern applications (see, e.g., [2]), which, de- 
pending on the level of detail of the knowledge to be represented, 
may feel happy with a high-level, qualitative language, or need to 
use a low-level, quantitative language. Qualitative languages suffer 

1 TCSPs stands for Temporal Constraint Satisfaction Problems, a well-known 
constraint-based temporal framework |6J. 

2 European Conference on Artificial Intelligence. 

3 The reviews are added to the actual paper, after the references, for potential 
people interested in objectivity of conferences' reviewing processes. 

4 Qualitative Spatial Reasoning. 



from what Forbus et al. 1 7 1 refer to as the poverty conjecture, but have 
the advantage of behaving computationally better. On the other hand, 
quantitative languages do not suffer from the poverty conjecture, but 
have a slow computatinal behaviour. Thus, such a KR system will 
feel happier when the knowledge at hand can be represented in a 
purely qualitative way, for it can then get rid of heavy numeric cal- 
culations, and restrict its computations to a manipulation of symbols, 
consisting, in the case of constraint-based languages in the style of 
the Region-Connection Calculus RCC-8 1 15 1, mainly in computing 
a closure under a composition table. 

An important question raised by the above discussion is clearly 
how to augment the chances of a qualitative/quantitative KR sys- 
tem to remain at the qualitative level. Consider, for instance, QSR 
constraint-based, RCC-8-like languages. Given the poverty conjec- 
ture, which corresponds to the fact that such a language can make 
only a finite number of distinctions, reflected by the number of its 
atomic relations, one way of answering the question could be to in- 
tegrate more than one QSR language within the same KR system. 
The knowledge at hand is then handled in a quantitative way only 
in the extreme case when it can be represented by none of the QSR 
languages which the system integrates. 

One way for a KR system, such as described above, to reason 
about its knowledge is to start with reasoning about the qualitative 
part of the knowledge, which decomposes, say, into n components, 
one for each of the QSR languages the system integrates. ForRCC-8- 
like languages, this can be done using a constraint propagation algo- 
rithm such as the one in Q. If in either of the n components, an 
inconsistency has been detected, then the whole knowledge has been 
detected to be inconsistent without the need of going into low-level 
details. If no inconsistency has been detected at the high, qualitative 
level, then the whole knowledge needs translation into the unifying 
quantitative language, and be processed in a purely quantitative way. 
But even when the high-level, qualitative computations fail to detect 
any inconsistency, they still potentially help the task of the low-level, 
purely quantitative computations. The situation can be compared to 
standard search algorithms in CSPs, where a local-consistency pre- 
processing is applied to the whole knowledge to potentially reduce 
the search space, and eventually detect the knowledge inconsistency, 
before the actual search for a solution starts. 

With the above considerations in mind, we consider the integra- 
tion of Frank's cone-shaped and projection-based calculi of cardi- 
nal direction relations 1 8 1, well-known in QSR. A complete decision 
procedure for the projection-based calculus is known from Ligozat's 



work [12 1. For the other calculus, based on a uniform 8-sector par- 
tition of the plane, making it more flexible and cognitively more 
plausible, no such procedure is known. For each of the two cal- 
culi, the region of the plane associated with each of the atomic re- 
lations is convex, and given by the intersection of two half-planes. 
As a consequence, each such relation can be equivalently written as 
a conjunction of linear inequalities on variables consisting of the co- 
ordinates of the relation's arguments. We consider a more general, 
qualitative/quantitative language, which, at the basic level, expresses 
convex constraints of the form r(x,y), where r is a cone-shaped or 
projection-based atomic relation of cardinal directions, or of the form 
(a, (3)(x, y), with a, (5 G [0, 2-k) and (/3 — a) G [0, 7r): the meaning 
of (a, /3)(x, y), in particular, is that point x belongs to the (convex) 
cone-shaped area rooted at y, and bounded by angles a and (3. We 
refer to such constraints as basic constraints: qualitative basic con- 
straint in the former case, and quantitative basic constraint in the lat- 
ter. A conjunction of basic constraints can be solved by first applying 
constraint propagation, based on a composition operation to be de- 
fined, which is basically the spatial counterpart of composition of 
two TCSP constraints 1 6 1. If the propagation detects no inconsisteny 
then the knowledge is translated into a system of linear inequalities, 
and solved with the well-known Simplex algorithm. The preprocess- 
ing of the qualitative component of the knowledge can be done with 
a constraint propagation algorithm such as the one in Q, and needs 
the composition tables of the cardinal direction calculi, which can be 
found in 1 8 1. 

The general form of a constraint is (si V • • • V s„)(x, y), which 
we also represent as {si, . . . , s n }(x, y), where Si(x, y), for all i G 
{1, . . . , n}, is a basic constraint, either qualitative or quantitative. 
The meaning of such a general constraint is that, either of the n ba- 
sic constraints is satisfied, i.e., si(x, y) V • ■ • V s n (x, y). A general 
constraint is qualitative if it is the disjunction of qualitative basic 
constraints of one type, cone-shaped or projection-based; it is quan- 
titative otherwise. The language can be looked at as the spatial coun- 
terpart of Dechter et al.'s TCSPs 1 6 1: the domain of a TCSP variable 
is IR, symbolising continuous time, whereas the domain of an SCSP 
variable is the cross product IR x IR, symbolising the continuous 
2-dimensional space. 

Due to space limitations, we restrict the presentation to the unify- 
ing TCSP-liks constraint language, will all the tools required to con- 
vince the reader that the work is the description of an implementable 
KR&R system, consisting in a search algorithm using path consis- 
tency as the filtering procedure during the search, and the Simplex 
algorithm as a completeness guarantee, at the leaves of the search 
space. 

2 Constraint satisfaction problems 

A constraint satisfaction problem (CSP) of order n consists of: 

1. a finite set of n variables, x\, . . . , x n ; 

2. a set U (called the universe of the problem); and 

3. a set of constraints on values from U which may be assigned to 
the variables. 

An m-ary constraint is of the form R(xi 1 , • • ■ , Xi m ), and asserts that 
the values Oj x , . . . , (Zj m assigned to the variables x% x , . . . , Xi m , re- 
spectively, are so that the m-tuple (a^ , . . . , fflj m ) belongs to the m- 
ary relation R (an m-ary relation over the universe U is any subset 
of U m ). An m-ary CSP is one of which the constraints are m-ary 
constraints. We will be concerned exclusively with binary CSPs. 



For any two binary relations R and S, R n S is the intersection of 
R and S, R U S is the union of R and S, R o S is the composition of 
R and S, and is the converse of R; these are defined as follows: 

Rf]S = {(a, b) : (a, b) £ R and (a, b) G S}, 

RUS = {(a, 6) : (a,b) G Rot(a,b) G S}, 

R o S — {(a, b) : for some c, (a, c) G R and (c, b) G S}, 

BT = {(a, b) : (b, a) G R}. 

Three special binary relations over a universe U are the empty re- 
lation which contains no pairs at all, the identity relation T\, — 
{(a, a) : a G U}, and the universal relation = U x U. 

2.1 Constraint matrices 

A binary constraint matrix of order n over U is an n x n-matrix, say 
£>, of binary relations over U verifying the following: 

(Vi < n)(Bii C J^j) (the diagonal property), 

(Vi,j < n)(Bij = (Bji)^") (the converse property). 

A binary CSP P of order n over a universe U can be associated with 
the following binary constraint matrix, denoted B p : 

1. Initialise all entries to the universal relation: (Vi, j < 

n)((B p )y ^-U) 

2. Initialise the diagonal elements to the identity relation: 
(Vi < n)((B p )u *~ T\j) 

3. For all pairs (an, Xj) of variables on which a constraint (xj, Xj) G 
i?is specified: {B p ) tJ <- (S P ) SJ n R, (B p ) jl <- ((S p ) lJ )". 

2.2 Strong ^-consistency, refinement 

Let P be a CSP of order n, V its set of variables and U its uni- 
verse. An instantiation of P is any n-tuple (01, 02, ... , o„) of U n , 
representing an assignment of a value to each variable. A consistent 
instantiation is an instantiation (01, 02, . . . , a„) which is a solution: 
(Vi, j < n)((ai, a,j) G (B p )ij). P is consistent if it has at least one 
solution; it is inconsistent otherwise. The consistency problem of P 
is the problem of verifying whether P is consistent. 

Let V' — {xi 1 , . . . , Xi . } be a subset of V. The sub-CSP of P gen- 
erated by V' , denoted Pi v t , is the CSP with V' as the set of variables, 
and whose constraint matrix is obtained by projecting the constraint 
matrix of P onto V: (Vk,l < j){(B P W' ) kl = (S p ) lfcl( ). P is k- 
consistent 1 9 10 :] (see also [4]) if for any subset V' of V containing 
k — 1 variables, and for any variable X £ V, every solution to P\yi 
can be extended to a solution to P\v'\j{x}- P is strongly fc-consistent 
if it is j-consistent, for all j < k. 

1-consistency, 2-consistency and 3-consistency correspond to 
node-consistency, arc-consistency and path-consistency, respectively 
1131 1141 . Strong n-consistency of P corresponds to what is called 
global consistency in |5|. Global consistency facilitates the impor- 
tant task of searching for a solution, which can be done, when the 
property is met, without backtracking 1 10 1. 

A refinement of P is a CSP P' with the same set of variables, and 
such that: (Vi,j)((B p ') !J C (B p )ij). 

3 A spatial counterpart of TCSPs: Spatial 
Constraint Satisfaction Problems (SCSPs) 

TCSPs (Temporal Constraint Satisfaction Problems) is a constraint- 
based framewrok well-known in Temporal Reasoning |6|. We pro- 
vide a spatial counterpart of TCSPs, which we refer to as SCSPs — 
Spatial Constraint Satisfaction Problems. The domain of an SCSP 



variable is the cross product IR x IR, which we look at as the set of 
points of the 2-dimensional space. As for a TCSP, an SCSP will have 
unary constraints and binary constraints, and unary constraints can 
be interpreted as special binary constraints by choosing an origin of 
the 2-dimensional space — space (0, 0). 

We first define some more terminology to be used in the rest 
of the paper. We make use of a Cartesian system of coordinates 
(0,x'x,y'y). The x-axis x'x is the origin of angles, and the an- 
ticlockwise orientation is the positive orientation for angles. Given 
that we use the set [0, 2-7r) as the universe of angles (measured in 
radians), if two angles a and (3 are so that a > /3, the interval 
('a, f3) 3 will represent the union Ifct, 2n) U [0, ft) 3 . Given a posi- 
tive real number a and a strictly positive integer n, we denote by 
a mod n the remainder of the integral division of a by n. Further- 
more, given any a, f3 G [0, 2tt), the difference (3 Q a will mea- 
sure the anticlockwise (angular) distance of (3 relative to a: i.e., 
P q a _ r P-a+2n mod 2j n . s i m ii ar iy ; the sum a © (3 of a and 

(3 is defined as a ® (3 = (2±£ mod 2)ir. 

Definition 1 (SCSP) An SCSP consists of (1) a finite number of 
variables ranging over the universe of points of the 2-dimensional 
space (henceforth 2D-points); and (2) SCSP constraints on the vari- 
ables. 

An SCSP constraint is either unary or binary, and either basic or dis- 
junctive. A basic constraint is (1) of the form e(x, y), e being equal- 
ity, or (2) of the general form If a, f3) J (x, y) (binary) or (fa, (3) J (x) 
(unary), with a, (3 G [0,2tt), {13 Q a) G [0,tt), i,j G {0,1}. '(°' 
and stand, respectively, for the left open bracket '(' and the left 
close bracket '['. Similarly, ') ' and stand, respectively, for the 
right open bracket ')' and the right close bracket ']'. A graphical il- 
lustration of a general basic constraint is provided in FigureQ 

A disjunctive constraint is of the form [Si V- ■ -V5„] (a;, y) (binary) 
or [Si V • ■ ■ V S n ](x) (unary), with Sk{x, y) and Sk{x), k = 1 . . . n, 
being basic constraints as described above: in the binary case, the 
meaning of such a disjunctive constraint is that, for some k = 1 . . . n, 
Sk(x, y) holds; similarly, in the unary case, the meaning is that, for 
some ft — 1 ... n, Sk(x) holds. A unary constraint R(x) may be 
seen as a special binary constraint if we consider an origin of the 
World (space (0, 0)), represented, say, by a variable xo: R(x) is then 
equivalent to R(x,xq). Unless explicitly stated otherwise, we as- 
sume, in the rest of the paper, that the constraints of an SCSP are 
all binary. 

An SCSP constraint R(x,y) is convex if, given an instantiation 
y — a of y, the set of points x satisfying R(x, a) is a convex subset 
of the plane. A universal SCSP constraint is an SCSP constraint of 
the form [0, 27r) (x,y): the knowledge consisting of such a constraint 
is equivalent to "no knowledge", i.e., any instantiation (a, b) of the 
pair (x, y) satisfies it. A universal constraint is also a convex con- 
straint. A convex SCSP is an SCSP of which all the constraints are 
convex. Given its similarity with an STP (Simple Temporal Problem) 
[6], we refer to a convex SCSP as an SSP (Simple Spatial Problem). 
An SCSP is basic if all its constraints are basic. We refer to a ba- 
sic SCSP as a BSP (Basic Spatial Problem). Note that a BSP may 
have pairs (x, y) of variables on which no constraint is specified (the 
implicit constraint on such pairs is then the universal relation [0, 7r), 
which we also refer to as ?). 

The standard path consistency procedure for binary CSPs is guided 
by three algebraic operations, the converse of a constraint, the com- 
position of two constraints, and the intersection of two constraints. 
These are defined below for SCSP basic constraints. The case of gen- 
eral (possibly disjunctive) constraints is obtainable from the case of 



y 

D' 




y' 

Figure 1. Graphical interpretation of the basic constraint ( l s, t) 3 (X, Y): 
Given Y, the set of points X satisfying the constraint ( l s, t) 3 (X, Y) is the 
cone-shaped area centred at Y, whose lower bound (open if i = 0, close 
otherwise) and upper bound (open if j = 0, close otherwise) are, respec- 
tively, the half-lines whose angular distances from the x-axis, with respect to 
anticlockwise orientation, are s and t. 

basic constraints. 

3.1 The converse of an SCSP basic constraint 

The converse of an SCSP relation R is the SCSP relation R~" such 
that, for all x, y, R(x,y) iff R~" (y , x) . We refer to the constraint 
R~"(y, x) as the converse of the constraint R(x, y). The converse of 
e(x, y) is clearly e(y, x). The converse of an SCSP basic constraint 
( l Q, (3) 3 (x, y) is the SCSP basic constraint ( l a ®tt,(3® Tv) 3 (y, x), 
which can be explained by the simple fact that, given any instan- 
tiation (x, y) = (a,b) of the pair (x,y) satisfying the constraint 
( I a, f3) J (x,y), the angle formed by the directed line (ba) with the 
x-axis is obtained by adding n to the angle formed by the directed 
line (ab) with the x-axis. 

3.2 The composition of two SCSP basic constraints 

Consider a point y of the plane, and an angle a in [0, ir). We de- 
note by l(y, a) the directed line through y forming angle a with the 
x-axis x'x. y and a partition the plane into five zones, which are 
the left open half-plane bounded by l(y,a), the half-line consist- 
ing of the points of l(y, a) coming before y (negative half-line), the 
point y itself, the half-line consisting of the points of l(y, a) com- 
ing after y (positive half-line), and the right open hal-plane bounded 
by l(y, a). We denote the five regions by lohp (y,a), nhl (y,a), 
pt-reg (y, a), phi (y, a), and rohp (y, a), respectively, and the set 
of all of them by REGIONS (y, a). Given a fixed angle a in [0, n), 
we can thus define a five-atom calculus CAC a of binary relations. 
The atoms are lohp a , nhl a , EQ, phl a and rohp a , defined as fol- 
lows, for all pairs (x, y) of 2D points: lohp a (x,y) iff x belongs 
to lohp (y, a), nhl a (x, y) iff x belongs to nhl (y, a), EQ(x, y) iff 
x = y, phl a (x, y) iff x belongs to phi (y, a), and rohp a (x, y) iff x 
belongs to rohp (y, a). We denote by ATOMS(a) the set of all five 
atoms. Clearly, lohp a (x, y) iff rohp a (y, x), nhl a (x, y) iffphl a (y, x), 
and EQ(x,y) iff EQ(y,x). In other words, lohp a and rohp a are 
each other's converses, and so are nhl a and phl a ; whereas EQ is 
its own converse. We consider now two fixed angles a and ft from 
[0, 7r) and compute the composition R\ o R 2 of Ri and i?2, with 
Ri G ATOMS(a) and R 2 G ATOMS((3). i?i o R 2 is the relation 
R — {(x, z) : for some y, Ri(x, y) and R 2 (y, z)}. Clearly, if Ri 



is EQ then Ri o R 2 = R 2 , and if R 2 = EQ then i?i o R 2 = R u 
We use the standard notation for (possibly) disjunctive relations. The 
other possibilities are presented in the (composition) table of Figure 
[2[Top), where: 
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lohpp 




phl fj 


rohpp 


lohp a 


ct 


eti 


ct 2 


? 


nhl a 


ctz 


Cti 


cts 


ct 6 


Phl a 


Ct\ 


Ct-j 


ct s 


Ctg 


rohp a 


? 


Ct a 


Ctb 


Ct c 



Ca,/jy s.t. 


Translation of ( l a, /3) J (x, y) 


a G [0,n),/3 G [0, tt) 


((lhp a yn(rhp g y)(x,y) 


a G [0,tt),/3 G [7T,27r) 


({lhp a y^(lhp g _^){x,y) 


a G [tt, 2tt),/3 G [7T,27t) 


({rhp a _ n yn(lhp g _„y)(x,y) 


a G [7r2,7r),/3 G [0,tt) 


((rh Pa _J>n(rhp fj y)(x,y) 



(3, whp a n rohp 3 
f3, rohp a n foftp^ 



Figure 2. (Top) Composition R o S, with t? atom of CA£ a and S atom 
of C4£/3. (Bottom) Translation of basic relation {'a, f3) 3 into R n S, with 
(possibly disjunctive) C4£ Q relation, and S (possibly disjunctive) CACp 
relation. 



• cto is lohp a if a = /3, ? otherwise; 

• eti is toAp^j if a > /?, ? otherwise; 

• ct2 is to/iPo, if a < /3, ? otherwise; 

• ct3 is ZoAp^j if a < /?, ? othrwise; 

• eti is /o/ip^ n rohpg if a > P, nhl a if a = (3, rohp a fl Zo/ip^ 
otherwise; 

• cts is /o/zp a n lohpg if a < f3, ? a if a = /3, rohp a n ro/ip^ 
otherwise; 

• ct6 is ro/ipfl if a > yS, ? otherwise; 

• ct7 is to/!p a n lohpg if a > f3, ? a if a = 
otherwise; 

• ctg is /oAp a (~l rohpg is a < /3, pW a if a 
otherwise; 

• ctg is rohpg if a < (3, phl a is a = /?, ? othrwise; 

• ct a is rohp a if a < [3,7 otherwise; 

• ctj, is rohp a if a > (3,1 otherwise; 

• ct c is rohp a if a — (3, ? otherwise; 

• ? = {(p, q) : p and g planar points} (i.e., ? is the universal binary 
relation on 2D points); 

• ? a = {(x,y) G? : x G %. a )} 
set of pairs (x, y) of 2D points s.t. x G /(y, a)) 

It follows from the above that, given a, f3 G [0, tt), the composition 
i?i o R 2 of i?i G ArOM5(a) and R 2 G Ar<9M5(/3) is a convex 
relation. 

It is now easy to derive the composition, R o S, of two SCSP basic 
constraints i? = { n a,f3) n and S = ( t2 -y,S) J2 . It is sufficient to 
know how to translate an SCSP basic relation R = ('a, f3) 3 into a 
conjunction R\ C\R 2 , where R\ isaC4£ a oraC4£ w _ a convex re- 
lation, and i?2 a C AC,3 or a C^4£^-^ convex relation: this is done in 
the table of Figure|2{Bottom), where the following notation is used. 
Given a G [0,n), we denote by lchp a (resp. rchp a ) the disjunctive 
relation {lohp a ,nhl a ,e,phl a } (resp. {nhl a ,e,phl a ,rohp a }). The 
constraint lchp a {x,y) (resp. rchp a {x,y)) means that x belongs to 
the Left (resp. Right) Close Half Plane bounded by l(y,a). Given 
a G [0, tt) and i G {0, 1}, the notation (lhp a ) 1 (resp. (rhp a ) 1 ) stands 
for lohp a (resp. rohp a ) if i = 0, and for /c/ip^ (resp. rchp a ) if i = 1. 
It is important to keep in mind, when reading the table of Figure 
|2jBottom), that a G [tv, 2tt) implies (a — tv) G [0, 7r). 



{nhl a , EQ, phi a } (i.e., the 



The composition Ro S of basic constraints R — ( ll a,/3)-' 1 and 
S = ( l2 7,5) J2 can thus be written as R o S = /(a) o 7(7) n 
/(a) o /(<5) n f((3) o /( 7 ) n /OS) o f(5), where /(*), for all x G 
{a, /3, 7,5}, is a C4£a; atom if x G [0,7r), and a CAC^-a, atom 
if 1 £ [7T, 27r). Given that, for all a, [3 G [0,7r), the composition 
i?i o R 2 of i?i G ATOMS(a) and 7? 2 G Ar<9MS(/3) is a convex 
relation, we infer that the composition of two SCSP basic constraint 
is an SCSP convex constraint. 

3.3 The intersection of two SCSP basic constraints 

Clearly, eHe = e; en('a,/3) J = eif 1 = j = 1; and en fa, /3) J = 
if z = or j = 0. 

Given a basic relation = ( l a,/3)-' and 7 G [0, 27r), 7 is an- 
ticlockwisely inside R (notation acwi(-y,R)) iff (1) 7 = and 
« = 1; (2) 7 = /3 and j = 1; or (3) 7 7^ a and 7 / (1 and 

/?ea=(79a) + (/?67)- 

It is now easy to derive the intersection, R fl S, of two 
basic constraints i? = and S = f 2 7,<5) J2 . If neither 

of acw^ajS), acwi(/3, S), acwi(-y,R) and acwi(<5, i?) holds, then 
R f] S = 0. Otherwise, the intersection is nonempty: R f~) S = 
0)* 7 ' If acwi(a, S) then 4> = a and « = ii, otherwise (j> — 7 and 
« = 12. If acwi(/3, S 1 ) then 9 = /3 and J ~ ji, otherwise 8 — S and 
J = J2- Clearly, if i? n 5 7^ then it is a basic constraint. 

The converse of an SCSP basic constraint is an SCSP basic con- 
straint. The composition of two SCSP basic constraints is either a 
basic constraint or the universal constraint. Finally, the intersection 
of two SCSP basic constraints is an SCSP basic constraint. Now, the 
only SCSP constraint that may (implicitly) appear in a BSP is, as 
already alluded to, the universal relation ?. Furthermore, the con- 
verse of ? is ?, ?n? =?, ?o? =?, and, for all basic relations R, 
Rf\1 =? n R = R and Ro? =? a R =?, This leads to the following 
theorem. 

Theorem 1 The class o/BSPs is closed under path consistency: ap- 
plying path consistency to fl BSP either detects inconsistency of the 
latter, or leads to a (path consistent) BSP. 

It remains, however, to be proven that path consistency terminates 
when applied to a BSP. Furthermore, if path consistency is to be used 
as the filtering method during the search for a path consistent BSP 
refinement of a general SCSP, then it should also be proven that path 
consistency terminates when applied to a general SCSP -it may be 
worth noting here that path consistency applied to a general TCSP 
1 6 1 may lead to what is known as the fragmentation problem 1 1 6 1 . We 
do this through the explanation of what we refer to as a "qualitative 
behaviour" of path consistency when applied to a general SCSP. 

3.4 Qualitative behaviour of path consistency 

Let P be a general SCSP and HOLES(P) the set of all 7 G [0, 2tt) 
such that there exists a constraint [Si V • ■ • V S n ](x, y) of P with, 
for some i G {1, . . . ,n}, Si of the form (^a.,f3y, and such that 
7 G {a, 13}. We also denote by HOLES + (P) the set HOLES(P) U 
{a G [0, tt) : (a + n) G HOLES(P)} U {a G [n, 2vr) : (a - n) e 
HOEES(P)}. Given a set A, we denote by \A\ the cardinality of 
A. Clearly \HOLES + (P)\ < 2 x \HOLES(P)\. The qualitative be- 
haviour comes from properties of the operations of converse, inter- 
section and composition when applied to SCSP basic constraints. The 
intersection R n S of two SCSP basic constraints R = ( ll a, f3) 31 
and S = C 2 7,5) J2 is of the form (*</>, 9 ) 3 , with both and 9 in 
{a, f3, 7, S}. The converse of an SCSP basic constraint ('a, /3) J is 





a = 


< a < | 


7T 

a = 2 


f < a < ty 


loh Pa (X,Y) 


yx > yv 


yx — yy > tga.(x x — Xy) 


yx < yy 


yx - yy > tg(-K - a).{xx - xy) 


lch Pa (X,Y) 


yx > yv 


yx - yy > tga.(x x - xy) 


yx < yy 


yx-yy > tg{n - a).(x x " x Y ) 


roh P JX,Y) 


yx < yy 


yx - yy < tga.(x x - x Y ) 


yx > yy 


yx — yy < tg(n - a).(x x - xy) 


rch Pa {X,Y) 


yx < yy 


yx - yy < tga.(x x - Xy) 


yx > yy 


yx — yy < tg(ir - a).(x x - Xy) 



Figure 3. Translation of an SCSP basic constraint into a conjunction of linear inequalities. 



('q © it, (3 © tt) 3 : but such an operation will not create new "holes", 
since if a basic constraint of the form { l a, f3) J (x,y) appears in P 
then we would have both a and /3 in HOLES(P), and both a © n and 
(3 ©7T in HOLES + (P). 

Path consistency using, as usual, a queue QUEUE where to put 
edges of P whose label has been updated, would thus, for each edge 
(i.e., pair of variables) (x,y), successfully update the label at most 
HOLES + (P) times. The number of edges is bounded by n 2 , n being 
the number of variables of P. Furthermore, when an edge is taken 
from Queue for propagation, 0(n) operations of converse, intersec- 
tion and composition are performed. This leads to the following the- 
orem stating termination, and providing a worst-case computational 
complexity, of path consistency applied to a general SCSP. 

Theorem 2 Applying path consistency to a general SCSP P with n 
variables terminates in 0(|HOLES(P) x n 3 ). 



high-level vision, such as, e.g., satellite-like surveillance of a geo- 
graphic area. We have provided all the required tools for the imple- 
mentation of the presented work; in particular, the algebraic opera- 
tions of converse, intersection and composition, which are needed by 
path consistency. An adaptation of a solution search algorithm, such 
as, e.g., the one in 1111 (see also which would use path con- 
sistency as the filtering procedure during the search, can be used to 
search for a path consistent BSP refinement of an input SCSP. But, 
because we know nothing about completeness of path consistency for 
BSPs, even when a path consistent BSP refinement exists, this does 
not say anything about consistency of the original SCSP. To make 
the search complete for SCSPs, we have proposed to augment it with 
the Simplex algorithm, by translating, whenever a leaf of the search 
space is successfully reached, the corresponding path consistent BSP 
into a conjunction of linear inequalities, which can be solved with 
the well-known Simplex algorithm 1 3 1. 



3.5 Translating an SCSP basic constraint into a 
conjunction of linear inequalities 

We now provide a translation of an SCSP basic constraint into (a 
conjunction of) linear inequalities. We will then be able to trans- 
late any BSP into a conjunction of linear inequalities, and solve it 
with the well-known Simplex algorithm (see, e.g., |3J). This will 
give a complete solution search algorithm for general SCSPs, us- 
ing path consistency at the internal nodes of the search space, as a 
filtering procedure, and the Simplex at the level of the leaves, as a 
completeness-guaranteeing procedure (the SCSP at the level of a leaf 
is a path-consistent BSP, but since we know nothing about complete- 
ness of path-consistency for BSPs, we need to translate into linear 
inequalities and solve with the Simplex). 

Given a point X of the plane, we denote by (xx,yx) its co- 
ordinates. The translation of e(X, Y) is obvious: xx — xy < 
A xy — xx < A yx — yy < A yy — yx < 0. For the translation 
of a general basic constraint ( l a, f3) 3 (X, Y), the results reported in 
the table of Figure[2lBottom) imply that all we need is to show how 
to represent with a linear inequality each of the following relations 
on points X and Y, where a G [0, it): lohp a (X, Y); lchp a (X, Y); 
rohp a (X, Y); and rchp a (X, Y). We split the study into four cases: 
a — 0, < a < ^, a = -|, 5 <a<7r. The result is given by 
Figure|5| where, given an angle a, tga denotes the tangent of a. We 
remind the reader that in a system of linear inequalities, there is a 
way of turning a strict inequality into a large one 1 3 1. 

4 Summary 

We have provided a TCSP-liks decidable constraint language for rea- 
soning about relative position of points of the 2-dimensional space. 
The language, SCSPs (Spatial Constraint Satisfaction Problems), 
subsumes two existing qualitative calculi of relations of cardinal di- 
rections 1 8 1, and is particularly suited for applications of large-scale 
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